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We study the probability distribution P of the sum of a large number of non-identically distributed 
random variables Um- Condensation of fluctuations, the phenomenon whereby one of such variables 
provides a macroscopic contribution to the global probability, is discussed and interpreted in analogy 
to phase-transitions in Statistical Mechanics. A general expression for P is derived, and its sensitivity 
to the details of the distribution of a single n„i is worked out. These general results are verified by 
the analytical and numerical solution of some specific examples. 
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I. INTRODUCTION 

Condensation is the phenomenon whereby a finite fraction of some quantity, e.g. a particle density, concentrates 
into a small region of phase-space, as in the paradigmatic example of a vapor transforming into a liquid when crossing 
a phase-transition. Condensation is observed in a number of different models, related to magnetic properties [I|, 
gravity mass transport and other issues 11®. Despite the prominent role played by molecular interactions in 

most cases, condensation can also be observed in non-interacting systems as, for instance, in a quantum Bose-Einstein 
condensate [2l| or in classical models such as the spherical model of a ferromagnet @ . In the above mentioned cases 
condensation occurs because the condensing quantity - the particle number in the former example or the spin variance 
in the latter - is conserved. This constraint acts like an effective interaction among the constituents bringing about 
the transition . Indeed, condensation does not occur in a non-interacting boson gas - as in the case of photons 

- which does not conserve the number of constituents. 

A different manifestation of condensation is observed when probability distributions of a fluctuating collective 
variable A/", such as the number of particles in a thermodynamic system, are considered. In this case, a fluctuation 
Af = N well above the typical value can be associated to a condensed configuration of the system ®-®|. This 
phenomenon, referred to as condensation of fluctuations, is not restricted to the particle number Af but was observed 
for quantities as diverse as energy, exchanged heats, particles currents etc... (22 - l3^ . It was shown ®, ®] that in 
some systems condensation of fluctuations may occur because, from the mathematical point of view, asking for a 
specific value Af = N constraints the system similarly to what a conservation law does. 

In this paper we study the probability distribution P of the sum Af oi M non-identically distributed random 
variables. We discuss how an interpretation can be provided, along the guidelines of Statistical Mechanics, in terms 
of a phase-transition between a normal phase with a vanishing order parameter and a condensed one with pi > 0. 
A general expression for the probability P is found and the radically different behavior of this quantity in the normal 
and in the condensed phase are discussed and illustrated by comparing the analytical and numerical solution of some 
specific models. In particular, in the condensed phase, the notable phenomenon of the giant response - a dramatic 
change of P as the statistical properties of even a single random variable is modified - is pointed out. 

This paper is organized as follows: In Sec. [Ill we introduce the statistical model that will be studied and set the 
notation. In Secs. nil and llVI its behavior is discussed when condensation does not occur and when it does, providing 
also an example by means of an analytically tractable case for identically distributed variables (Secs. IIII Al and llV A[ 
respectively). In Sec. |V]the case of non identically distributed variables is addressed and the phenomenon of the giant 
response is discussed (Sec. IV Al) . Some examples are considered in Sec. IV Bl Finally, In Sec. I VII we briefly summarize 
and conclude the paper. 


II. THE STATISTICAL MODEL 

In order to set the stage, let us consider the independent variables Um {m = I,M) subject to a probability 
Pm{nm] K), where AT is a set of parameters. This probabilistic setup is suited to describe at a simple level a variety 
of systems ranging from physics to chemistry, biology and social sciences. For instance, one can consider M receptors 
where Um ligand particles, like those of a pollutant, can be adsorbed, or Um electrons populating M atomic levels. 
One can also think of M individuals, or agents, collecting Um resources with a certain probability Pm- In the former 
examples the temperature can be one of the control parameters K but, in general, others can be present. 
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We fix the language by speaking of M receptors hosting a total number 

M 

Af ='^ Tim ( 1 ) 

m—1 


of particles, with an average value {Af) = '"^here (rim) = '^,^npm{n',K). For ease of notation the 

dependence on K will be often dropped, and M will be considered large. 

We are interested in the probability to observe a total number Af = N oi particles 


P{N, M) 


Vi{ni)p2{n2)-■ ■pMinM)SM,N = 

ni,n2,...,nM 


1 

27ri 




( 2 ) 


where, for discrete variables, we used the representation 5f,f,N = ^ § dz z A^+i)^ with 


Q{z,M) 


■ M 

_m—l rim 


1 

M 


( 3 ) 


and p = ~ ^ is the particles density. 

As explained in [3, 1^ . l23j | the probability distribution P of the fluctuations of the particle number corresponds 
also to the partition function of a dual model where the number of particles is conserved. For example, with some 
particular choices of the microscopic probabilities p that will be considered below, such dual model corresponds to 
specific instances of the so called urn model (or balls in boxes model) or of a zero-range process. This duality, which 
to the best of our knowledge was never discussed in connection to the above mentioned models, allows us to borrow 
a number of well established results in this research areas to illustrate the behavior of P - whose properties are here 
discussed in a rather large generality - in some exemplifying cases. 

The following relation holds [ 1 , @ 


N 

P{N,M) = E TT{n,N,M) 

n—0 


( 4 ) 


where 


Tr{n,N,M) = P{N — n,M — l)pM{n) (5) 

is the probability that, when the M-th receptor is added to the previous M — 1, n particles are stored in it. eq. ® 
is a recurrence relation allowing one to determine the probability distribution of M variables once the one for M — 1 
is known. 

Let us discuss the basic mechanisms whereby the recurrence (|3]) works to build up the full probability P{N,M). 
Denoting with N = Mp the value of N where P is maximum, the largest P{N — n, M — 1) in eq. ([5]) is the one with 
n = N — N. Notice that this term is present in eq. ® only if > TV. 


III. GAS PHASE 

According to the previous discussion, for TV < TV, tt has a maximum at a value n = rig which is microscopic and 
does not scale with M (see the inset of figlDfor a specific example with rig = 0, to be discussed below). This is so 
because the quantities P(TV — n,M — 1) in eq. (I5|) lower with n (i.e. moving away from the maximum) and the same 
is true for pM(ri), for sufficiently large n > rig (being pMin) normalized). If the Pm’s are monotonously decreasing, it 
is rig = 0. 

A similar setting, with tt peaked at a microscopic rig, is found also for N > N, when the largest probability 
P{N,M — 1) is contained in the sum on the r.h.s. of eq. (|T|), but its maximum is tamed by the microscopic 
probabilities, i.e. limM->.oo P{N, M — 1)pm(TV — TV) = 0. 

The situation with tt peaked in n = rig is physically intuitive: It expresses the fact that, when M is large, the 
occupancy tt of the new receptor (the TH-th) is microscopic. We will denote this situation, with a uniformly small 
occupation, the normal (or gas) phase. 
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A. An example 


Let us illustrate these behaviors by considering a specific example with power-law distribution 

Pm{n) = + ( 6 ) 

where Km > 1 and c is a normalization. 

We start with the simplest case where Km = K does not depend on m [s^- In the inset of figlU tt is plotted for 
K = ^j2 and different choices of N. Here one observes a sharp peak in n = = 0, as expected. 



FIG. 1: The quantity n{n, N, M) is plotted against n for the model with probabilities ((6]) with a uniform Km = K = 3 > Kc, 
for M = 246 and the values of N indicated in the key. The location of Ugi is indicated with a bold vertical segment. In the 
inset the same plot is made for Km — K = 3/2 < Kc- 


In the gas phase, P can be determined in the large-M limit by a steepest-descent evaluation [3, H, 0, [TtI - I^ . [H, 

[s^, [ 3 ^ of the integral on the r.h.s. of eq. leading to 

~ (7) 

with an M-independent rate-function 

R{p) =-\nQ[z*{p)]+plnz*{p), (8) 

where z*{p) is given by the saddle-point equation 


z*{p) 


Q'[z*ip)] 

Q[z*{p)] 


= P- 


(9) 
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For the specific example above, this equation can be cast as 


LiK-i{z*) 

LiKiz*) 


p + 1, 


( 10 ) 


where Li^iz) is the polylogari thm (Jonquiere’s function). It admits a solution with z < 1 for any finite value of p if 
K <K^ = 2 1 P is then expressed by Eqs. (I7I8D with 

Q{z*)=Lik{z*). (11) 

The rate-function R{p) obtained in this way is plotted with a black heavy dashed line in figlU In the same picture 
the quantity 

n{p,M) = -^\nP{Mp,M) (12) 

obtained from eq. ([2|) by exact enumeration, is shown for different choices of M. One observes that TZ approaches 
the asymptotic M-independent form R{p) as M is increased. Notice that the convergence is faster at small densities. 
It must be recalled that, for K < Kc, the average number (N) of particles is not finite, meaning that in the large-M 
limit fluctuations with large N are very likely, as it is reflected by the vanishing of R{p) at large densities. However, 
for finite M such large values of N cannot be sustained and TZ, after reaching a minimum, raises again increasing p, 
thus determining the existence of a most probable value of the fluctuations. The position of such value is pushed to 
larger densities by increasing M, providing in this way a gradual convergence, from smaller to larger values of p, of 
TZ towards R. 



P 


FIG. 2: 7Z(p, M) is plotted against p for for the model with probabilities (O with a uniform = K — 3/2 < Kc, for different 
values of M in the key. The heavy dashed line is the asymptotic (for M —>■ oo) expression R{p) obtained from Eqs. (I8I10I11II . 
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IV. CONDENSED PHASE 

A radically different situation occurs for N > N, when P{N, M — 1) grows fast enough to give liniM->-oo P{N, M — 
— A^) 0- For the choice ([5]) with = A', this happens when K > K^- 

In this case the sum in eq. (U) does not only take contributions around the microscopic value Ug, since tt can be 
non-negligible up to a certain n = ni oi order M . This is because tt is triggered by the maximum of P{N — n, M — 1) 
at a value n = N — N which, for given p, is of order M itself. Usually P is sharply peaked around the maximum as 
to give ne ^ N — N for N sufficiently larger than N. FigH] shows this for the model ([6|) with A" = 3: tt develops a 
second peak a,t n = m ~ N — N and only for n > it drops down to negligible values. The properties of this second 
peak can be described using extreme-value statistics 

The physical interpretation is that, when a value of N larger than the typical one N is attempted, the occupancy of 
the new receptor M can be either microscopic or macroscopic. Then, together with the uniformly scarcely populated 
gas phase, a liquid, or condensed phase coexists characterized by a single Um hosting a finite fraction of the N particles. 

From the previous considerations a close resemblance emerges with the problem of a gas-liquid transition, with N 
being a control parameter playing the role of the volume, pi = ni/M an order-parameter and H ~ — InTr an energetic 
landscape. Notice that the amount of condensed vs normal fluctuations, described by pi, depends on N: The fraction 
of condensate is absent a,t N = N, and increases with N. 

When condensation occurs eq. (fTOl) has no solutions. Therefore the steepest-descent evaluation of the integral in 
(I2|) cannot be carried over straightforwardly as done for the gas phase in Sec. IIIII In [s^ an upgraded saddle point 
technique based on a density functional approach was shown for a case with continuous variables related to the example 
®. Another way of proceeding - still resorting to the saddle point technique - is illustrated in the Appendix. Here 
we prefer to determine the form of P in a different way, which is discussed now. Starting for simplicity from equally 
distributed variables, Pm{n) = p(n), we assume that particles condensed in a single receptor contribute a fraction 
C{N, M) = (1 — a)P of the global probability while the others, scattered over the remaining locations, provide the 
remaining part G{N, M) = aP. The parameter a depends on N in such a way that a = 1 (there is no condensate) at 
N < N, while a —)■ 0 (all the particles are condensed) for N ^ oo. Casting eq. ([4]) as 

rigt N 

P{N,M)^'^TT{n,N,M)+ 7r{n,N,M)=G{N,M) + C{N,M), (13) 

n—0 n—rigi 

where Ugi is the value of n where tt is minimum (see figHI), allows one to identify G and G as the first and second 
sum on the r.h.s., respectively. 

For N ':$> N the peak around ni becomes sharper and a gaussian approximation for the evaluation of the second 
term gives 


G{N, M) = (1 - a)P{N, M) ~ aP{Ni, M - l)p{N - Ni), (14) 

where Ni = N — ni N and cr = •\/27r[— (9^ hnr/ dn^)n=nt\~^ ■ Assuming that in the large-N limit aP{Ni, M) has 
only a weak dependence on N eq. (jI4|) has an approximate solution with 

aPiNi,M -l)c^il-a)bM, (15) 

and 

P{N,M) ~bMp{N -N). (16) 

where we have confused Ni with N for large N as expressed below eq. m- Eqs. (1161) is a general expression for the 
probability in the condensed phase for identically distributed variables. A straightforward generalization to the case 
of non-identically distributed variables will be discussed below. 

The quantity 5 m depends on the structure of the microscopic probabilities. If the p’s do not depend on M, we 
can infer 5 m by observing that ~ number of particles in the condensed phase - is of order M. 

Recalling Eqs. (I14I15I16|) one can argue that 


5m M. (17) 

The results (I16I17|) have a very transparent physical meaning [^, : when condensation occurs a single receptor 

hosts a number N — N oi particles with a probability p{N — N). The factor M is the number of ways to choose such 
receptor out of M, which is true when the receptors are identical (indeed we will show in Sec. 0that eq. (flT)) can be 
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violated for non-identically distributed variables). Notice that our derivation is of a general character and does not 
rely on any specific form of the microscopic probabilities p{n). 

We emphasize the crucial role played by the 5-function in eq. ([2|) which, as mentioned in the introduction, effectively 
constraints the total particle number thus invalidating the central limit theorem which would otherwise apply for the 
problem at hand, making the condensation phenomenon possible. 

Before moving to the more general case of a non-identical distribution of the Pm’s, we illustrate all the above with 
an example. 


A. An example 


Let us consider again the distribution ([6]) with K > Ke¬ 
rn exists only for p < p defined by 


It is easy to show that the saddle point solution to eq. 


LiK-i{l) 


(18) 


[p ~ 0.368433 for K = 3). For p < p only the normal phase exists, the steepest descent evaluation of the integral in 
eq. m is appropriate, and one arrives at eqs. (I8I10I) . The rate function obtained in this way is shown in the upper 
panel of figE] together with the behavior of TZ{p,M) [eq. (fT^ ] which, for p <'p, approaches R{p) for large M. 

As already discussed, for values of the density larger than p a straightforward steepest-descent evaluation of the 
integral in ([2]) breaks down. In this case P is not exponentially small in M, as required by eq. (ED, as it can clearly 
be understood observing in fig. [3] that the dependence on M does not cancels and TZ keeps decreasing to zero for any 
value of M. In this region condensation occurs and, instead of eqs. ()8I10|1 . the solution ()16I17I1 applies. In order to see 
this, in the inset of figElwe plot M~^P{N, M), since according to eq. (I16I17P this quantity ought to be independent 
of M and proportional to p{N — N). As expected, for p ^ the form (I16I17|) describes the probability with great 
accuracy for large N. Notice that in the condensed region the convergence to the asymptotic form [eq. (I16I17D ] is 
much faster than the one in the gas phase [to eqs. EHHD], being achieved already for M < 100, a value for which TZ 
is still quite different from R in the region p <])■ 


V. NON IDENTICALLY DISTRIBUTED VARIABLES 

Now we turn to study the phenomenon of condensation when the microscopic variables are not identically dis¬ 
tributed. Specific issues of this and related problems have been addressed in |3. fl^. I36l - l4^ . Here we are interest in the 
derivation of a general form for P, generalizing eq. (US, and to discuss the related phenomenon of the giant response 
on broad grounds. 

When the variables are non-identically distributed one can argue that condensation occurs on the most favorable 
receptor @,113, namely the one with the larger Pmini). In the example (jG]), it is the one with the smaller Km. 
Denoting m this term [i.e. Pm{ng) > Pmin-g) Vm ^ to ], recalling the physical meaning of tt in eq. ED, it is clear that 
a structure like the one in Fig. El with a sharp peak around will be present if the recently added receptor is the 
one where condensation occurs, namely if pm = Pm in eq- ED- Then, in order to proceed as in Sec. IlYl we define 

-K{n,N,M) = P{N-n,M-l)pj^{n), (19) 

which amounts only to the choice of a particular labeling of the receptors. Proceeding as in Sec. IIVI one obtains the 
following equation 


C(V, M) = (I - a)P(V, M) - <jP{Ng, M - l)pm{N - Ne), 
instead of eq. El, thus arriving at 


( 20 ) 


P{N,M)c^bMPm{N-N) (21) 

in place of eq. El- This form of the probability generalize eqs. (I16I17I) to the case of non-identically distributed 
variables. Notice that no assumptions on the form of the microscopic probabilities Pmin) has been made also in this 
case and, therefore, eq. m is expected to hold quite generally. Together with eqs. (I16I17I) . this equation represents 
the main result of this paper. Notice that the dependence on M of 6 m can be very different from the one EZD holding 
for identically distributed variables. Indeed, when the microscopic probabilities depend on M, the number of particles 
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in the condensed state may not be simply proportional to M, since the M-th receptors can promote condensation 
differently from the previous ones. An example showing this will be shown in Sec. IV Bl A straightforward consequence 
of eq. (l2T]l is the phenomenon of the extreme sensitivity of the global probability P to specific details of the microscopic 
ones Prm that we discuss below. 


A. Giant response 

It must be stressed that the solutions (I16I21I) are totally different from the one ([71), in particular concerning the 
dependence on M. Indeed, while (O is exponentially small for large M and transforms into a 5{N — N), eq. (I16I17I) 
shows that in the presence of condensation the dependence can be as weak as linear in M, signaling the occurrence 
of anomalously large fluctuations. 

Related to that, an extreme sensitivity of the macroscopic probability P to the details of the microscopic ones Pm 
arises. In fact, eq. (|2T]) clearly shows that the distribution of the single variable m fully determines the global quantity 
P. Introducing a susceptibility x ~ the shift of the macroscopic probability due to the variation of a microscopic one 
- from eq. (EU one has 


X{N, m) 


lim 

M —¥00 


AP{N, M) 

Apra{N) 


0, gas 

„ condens. 
PMdm,rn, 


( 22 ) 


where Bm is a constant. This shows that in the gas phase a shift of one (or even of a finite number) of the microscopic 
probabilities cannot alter the global behavior of P, since this is determined by the synergic contribution of a number 
M —> oo of variables, the situation is profoundly different when condensation occurs. In this case P is fully determined 
by the statistical properties of the most favorable receptor m. Therefore P is independent of the form of all the other 
AI —1 variables, whereas a macroscopic effect can be determined by altering the statistical properties of the single 
As we will show by means of some examples in Sec. IV Bl this may have dramatic effects on the form of the probability 
P in the condensed phase. 

This anomalous susceptibility is reminiscent of the large response induced by gapless modes [i^, such as massless 
Goldstone modes in systems with a spontaneously broken continuous symmetry. Actually, starting from equally 
distributed variables the symmetry between the receptors is broken by changing the properties of one of them. 
However the phenomenon of the giant susceptibility discussed here is more general since it occurs also when the 
modification of the probability of a single variable occurs in a set of (already) non-identically distributed ones, as will 
be illustrated by the second example of Sec. IV Bl (fig. H]). To the best of our knowledge, this remarkable property of 
the susceptibility (1771) was never pointed out before. 


B. Examples 

The occurrence of condensation in the case of non-identically distributed variables and the giant response phe¬ 
nomenon can be illustrated using again the probabilities (jH]) with m-dependent K^s. The simplest non-trivial choice 
is when all the K^s are equal except one, namely Km = K, Vm > 1 while Ki can be different from K. 

In the lower panel of figI3]we compare P for the three cases i) Ki = K = 3, ii) Ki = K = 6 and hi) Ki = 3, AT = 6. 
One sees that the curves relative to the choices ii) and iii) coincide for p <]} (i.e. up to the maximum of P). This 
is because in this region there is no condensation and the macroscopic probability is insensitive to a single pm, eq. 
(ITT)) . However, for p > p the two curves become totally different and instead case iii) behaves as i), apart from a 
vertical displacement due to the different value of the constant 6 m in eq. m- This shows that a single variable 
cannot influence the collective behavior unless it is the one where condensation occurs, in which case a giant response 
is observed. 

The examples considered insofar where based on power-law probabilities (jS]). However, the features above are more 
general and not only restricted to this case. We show this by considering the exponential form 

Pm(n) = c exp [-A"„ • n'"™], (23) 

where Km = P iP ^nd c are constants), and Km an M-independent exponent. This case is interesting also 

because the microscopic probabilities p do not depend only on m, but also on M. In this case the scaling dnD, which 
was expected quite generally for M-independent p’s, can - in principle - be spoiled, and a general form of 6 m is not 
available. 












In order to illustrate the behavior of P with the exponentially distributed microscopic probabilities we have evaluated 
it for different choices of the parameters entering eq. (USD- Starting with a uniform exponent = k = 1, setting 
/3 = 1 and a = 2, the upper panel of hgH] shows a pattern of behavior similar to the case ([HI) with K > Kc- for 
N < N, TZ approaches the form © with a rate function given by Eqs. (E©), whereas for TV iV the determination 
_ (tv-77) 

(ED holds (with m = 1), implying P{N,M) ^ e , as shown in the inset. The data collapse is obtained by 

plotting M‘^P{N, M) against (TV — N)/M‘^, implying that 6 m = bM~‘^. Notice that the approach to the asymptotic 
form is much faster than for the fat-tailed probabilities ©, since already for M ~ 100 one has a good representation 
of the large-TVf form in the range of densities considered, at variance with what observed in Fig. [3] 

The phenomenon of the giant susceptibility is illustrated by comparing the case above with the one where we change 
the distribution of ni as to have ki = 0.95 and all the remaining ones are left untouched (k^, = n, Vk > 1). The lower 
panel of figH] shows that, while in the normal phase TV < TV this does not alter P (a residual difference between the 
two curves is due to the finite value of TkT), a dramatic change is produced in the condensed region TV > TV because, 
since ki < k, the statistical properties of the condensing variable have been changed. 


VI. CONCLUSIONS 

In this paper we have discussed the general problem of evaluating the probability distribution P{N, M) of the sum 
TV of a large number M of micro-variables, not necessarily identically distributed. 

We have done this by means of the recurrence relation (©, which provides an analogy with a thermodynamic 
system where a condensation transition occurs and the identification of an order parameter pi. Eq. I© allows also the 
derivation of a rather general expression for P [eq. (1211) ] which is valid, when condensation occurs, for finite values of 
M . From this expression, computing the susceptibility (1221) . the extreme sensitivity of P to the distribution of even 
a single variable was explicitly shown. 

These properties of the probability P have been discussed by means of specific examples amenable of analytical 
and numerical computations, including identically and differently distributed variables, with or without fat-tails and 
also in the case of a specific dependence of the microscopic probabilities p on the number M of micro-variables. 

The noteworthy features discussed in this paper are associated to the existence of a condensation phenomenon 
and, therefore, they are not expected to be only relevant to the large deviations of A/", but also to those of different 
macrovariables, and to apply to a large class of problems in Physics and other areas, making the issue considered in 
this paper a broad and general research topic. 
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VII. APPENDIX: SADDLE-POINT EVALUATION OF P IN THE CONDENSED PHASE 

In the condensed phase the symmetry between the receptors is broken since, as discussed in Sec. IIVI a single 
variable provides a contribution comparable to all the remaining ones. Let us, without loss of generality, indicate this 
variable as being the first, namely ni. In view of that, we re-write the first line of eq. © as follows 

P{N,M) = ^p(ni)H(ni,TV) (24) 

ni 

with 

fl{ni,N,M)= ^ P 2 in 2 )---PMinM)SAf^,N-m, (25) 

n2,n3,...,nAf 

where A/) = X]m=2 ^ solution can be found by making the ansatz that in the condensed phase the argument of 

the sum in eq. (1241) is sharply peaked around a certain value ni, with a certain width b, so that it can be evaluated as 


P{N, M) = Mbp{ni)n{ne, N, PI) 


(26) 
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where b = Y^27r[—(9^ and the factor M in front of the r.h.s. of eq. (|26l) is due to the M possible 

ways of choosing the variable denoted by ni among M. Using eq. (IMl) as a starting point, instead of eq. , one 
arrives at 


P(N,M) =&Mp(nU— 

27ri J 


(27) 


for large M, where now 


Q{z,M) 


■ M 

_m—2 rim 


1 

M-l 


(28) 


and Pi = ni/M is the condensed particles density. The steepest descend evaluation of the integral leads to the saddle 
point equation 


^ Q ^ (p) 

^ (P) nr *( M = P' 


Pi- 


(29) 


In the condensed phase there is always a solution with z = 1 and pi = p — , and the evaluation of the integral 

in eq. gives 


P{N,M) = bMP{ni) (30) 

with 6 m '^bM /(2Tri) = bM/{2'iTi) (in the last passage we have used <5(1) = 1 because of the normalization 

of the p’s). Recalling that ni N — N (see Sec. IIIip one recovers the result (fT6)l that was obtained in a different way 
- by using the recurrency relation (|31) - in Sec. uni 
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p 



FIG. 3: Upper panel: The function TZ{p, M) is plotted against p for the model with probabilities (jb]) with a uniform Km = 
3 > Kc, for the values of M in the key. The heavy-dotted black line is this quantity for M —>■ oo, which coincide with R{p) 
obtained from Eqs. (I8I10II for p < ]} and is identically zero for p > 'p. In the inset M~^P is plotted against N — N. The 
heavy dotted turquoise line is the law {N — N)~^, i.e. eas. (ll6ll7ll . Lower panel: P is plotted for M — 333 and the three 
different choices (see text) i) Km = K = 3, Vm, continuous brown, ii) Km = K — 6, Vm, green dashed with a square and iii) 
Ki = 3, Km — K — 6, Vm > 1, dotted maroon with a circle. 
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FIG. 4: Upper panel: The function TZ{p, M) is plotted against p for the model with exponential probabilities (see text) with 






















